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ABSTRACT 

The Sunyaev-Zel'dovich (SZ) effect is the inverse Compton-scattering of cosmic microwave background 
(CMB) photons by hot electrons in the intervening gas throughout the universe. The effect has a distinct 
spectral signature that allows its separation from other signals in multifrequency CMB datasets. Using CMB 
anisotropies measured at three frequencies by the BOOMERANG 2003 flight we constrain SZ fluctuations in 
the 10 arcmin to 1 deg angular range. Propagating errors and potential systematic effects through simulations, 
we obtain an overall upper limit of 15.3 /j,K (2a) for rms SZ fluctuations in a broad bin between multipoles of 
250 and 1200 at the Ray leigh- Jeans (RJ) end of the spectrum. The resulting upper limit on the local universe 
normalization of the density perturbations with BOOMERANG SZ data alone is cr| z < 1.14 at the 95% con- 
fidence level. When combined with other CMB anisotropy and SZ measurements, we find cr| z < 0.92 (95% 
c.l.). 

Subject headings: cosmic microwave background — cosmological parameters — cosmology: observations — 
large-scale structure of universe 



1. INTRODUCTION 

The Sunyaev - Zel'dovich (SZ) effect dSunvaev et al. I 
1 19721) is scattering of cosmic microwave background (CMB) 
photons by electrons in the intervening gas throughout the 
universe. The SZ effect is generally subdivided into two 
subcomponents: the kinetic effect due to the bulk mo- 
tion and the thermal effect due to energy transfer from 
hot electrons. The latter is expected to be the largest 
modification to the background temperature during pho- 
ton transit from the last scattering surface. We will con- 
centrate on it hereafter and refer to it as the SZ sig- 
nal. The strongest SZ perturbations to the CMB temper- 
ature are dominated by scatteri ng of photons via hot elec- 
trons in massive galaxy clusters dKomatsu & Kitavamall 19991: 
ISpringel et al~l 120011: ICooravl 120001: iMolnar & Birkinshawl 
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12000: ISeliak et al. Il200lt ISadeh & RephaellloM .faiyhgre 
could also be an SZ signal from reionization (lOh et al. I 
120031) . The integrated SZ angular power spectrum is now a 
known probe of the amplitude of densi ty perturbations, a% 
dKomatsu et al. Il2002t Bond et al. 112005b . 

The SZ effect has been clearly imaged toward individual 
galaxy clusters dGreggo et al. Il2000l ICarlstrom et al. 1 120021: 
iJones et al. I Il993h~ and has been used for a variety of ap- 
plications, including a measurement of the CMB tempera- 
ture a t the redshifts of Coma and A2163 (Battistelli et al. 
|2002|) . However, the amplitude of the SZ power spec- 
trum at arcminute angular scales, generated from unresolved 
galaxy clusters, is still not well established with differ- 
enc es at the 2a level from a variety of detections and lim- 
its dSievers et al. ll2009HReichardt et al. Il2008l:lpawson et al. 
20061: ISharp et al. 1120091: iFriedman et al. ll2009HSayers et al. 
20081) . Existing SZ anisotropy measurements are restricted to 
observations with a narrow frequency coverage and to small 
areas on the sky. The differences could be a combination of 
foreground co ntamination an d large non-Gaussian variance of 
the SZ signal dCooray||2001l) . Also, no constraints on the SZ 
signal exist at tens of arcminute scales where primary CMB 
fluctuations dominate. 

A clean separation of the SZ anisotropies from primordial 
CMB is possible due to the fact that the SZ signal has a 
distin ct frequency spectru m from the 2.7K blackbody spec- 
trum dCooray et al. Il200"ol) . The spectral difference arises as 
inverse-Compton scattering leads to, on average, a net energy 
gain for the CMB photons and the scattered photons move 
from the low frequency Raylei gh-Jeans (RJ) tail to high fre- 
quencies dSunvaev et al. Il972h . The SZ sky is colder than the 
CMB at low frequencies and hotter than the CMB at high fre- 
quencies with no difference at about a frequency of 217 GHz. 
A potential detection of the large-scale structure SZ fluctua- 
tions is then aided by observations across the SZ null from the 
negative side to the positive side. 

A dataset of the form needed for a study of the large- 
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covariance of data multipole moments C : 



*lm u lm 



subject to 



0-0 0-1 0-2 1-1 1-2 2-2 



FIG. 1. — Elements of the matrix Cy (Eq. 1) for the data (blue triangle) 
and simulations in (/^K) 2 RJ units in the first multipole bin (250 <l< 450). 
The six independent points of the Cn matrix correspond to correlation be- 
tween frequencies labeled 0, 1, and 2 for 145 GHz, 245 GHz, and 345 GHz 
respectively. We show two sets of distributions for simulations. The smaller 
distribution takes into account only instrumental noise and cosmic variance. 
The larger distribution includes also primary CMB and foregrounds by tak- 
ing the rms of three different amplitudes for the components (see text). In the 
case where only one error bar is visible, the rms from foreground is smaller 
than the instrumental noise and cosmic variance. The inserted plot shows the 
distribution of Coo with the data value inserted as a vertical blue line. 

scale structure SZ effect is provided by the 20 03 flight of 
the b alloon-borne BOOMERANG experiment ( Masi et al. I 
120061) . This instrument derives directly from the BOOM- 
ERANG pay load that was flown in 1998 and resulted in first 
high signal-to-nois e maps of the CMB anis otropy with sub- 
horizon resolution dde Bernardis et al. 120 00). The instrument 
was launched by NASA on 2003 January 6 from Williams 
Field near McMurdo Station, in Antarctica. The flight lasted 
a total of 311 hours until 2003 January 21 and 119 hours of 
this observing period were devoted to scanning a deep survey 
region. The remaining time was spent on scanning a larger 
shallow survey and a section of the Galactic plane. Here, we 
concentrate on a search for SZ effect in the central deep field 
over 100 deg 2 with the highest signal-to-noise ratio. 

In this Letter, we report the first statistical limits of the SZ 
signal at sub-degree angular scales at these wavelengths. The 
discussion is organized as follows: in Section|2]we detail our 
approach to extract the SZ signal from the multifrequency 
BOOMERANG dataset; in Section [3] we detail our simula- 
tions used to estimate statistical and systematic uncertainties; 
and in Section|4] we present our results. 

2. CMB AND FOREGROUND REMOVAL 

To separate the SZ signal from all other sources 
of anisotropies, we adopted a technique well-known in 
the literature for removing foregrounds f rom the CMB 
aniso t ropies dTegmark & Efstathioul 119961 : iTegmark et al. I 
l200l lAmbTard et al. H2007T In our case, instead of recov- 
ering the primordial CMB signal, we recover SZ fluctuations 
by minimizing the covariance relative to the SZ frequency de- 
pendence, and treat primordial fluctuations as another source 
of noise. In the remainder of this Letter, we will refer to "fore- 
grounds" as all the additional emissions (CMB, Galactic dust, 
far-IR sources or FIRB, radio point sources) to the SZ effect. 

The power spectrum of the SZ can be obtained from a 
weighted mea n of the power spect ra at different frequencies: 
Cf z = w r Cw dCoorav et al. 1120001) . The weights w\ at each 
frequency i and multipole I can be obtained by minimizing the 



the constraint that SZ estimation is unbiased (j^-w^j) = 1). 

In the case of BOOMERANG maps, each frequency chan- 
nel consists of several detectors. We distinguish individual 
detectors with the indices (a, (3), while are indices for 
the frequency channels. In order to minimize instrumental 
noise more aggressively, we compute the covariance matrix of 
the signal by averaging all the combinations of cross-spectra 
between different detectors and ignoring the auto-spectra of 
the same detector. The contribution of the correlated noise 
between two different detectors is taken into account in the 
simulations and removed as part of a residual contribution to 
the SZ signal. 

We construct the binned covariance matrix in multipole I 
bin b as 



I i,a j,f) 



U3 



withal/3, if i= j , (1) 



where s(z/,) is the SZ frequency dependence at each of the 
BOOMERANG frequency bands relative to CMB with s(v) = 
2-(^/2)coth(jc/2), x = Iw/kTcMB ~ ^/56.8GHz, and b]' 01 is 
the measured beam window function for the detector a in 
channel i. Note that with the definition above, in the RJ limit 
s{v) — > 1 so that Cf z O, v 1 ) = s{v)s{v')Cf z where Cf z is the 
SZ anisotropy power spectrum in the RJ limit. The covariance 
matrix Cy is required to be invertible and positive definite. We 
numerically check this both in data and simulations. In Fig- 
ure Q] we show the covariance matrix from data and compare 
it to simulations described below. 

Using the data covariance matrix, the optimal weights for 
the SZ reconstruction are 



C-'e 
e r C-'e 



(2) 



where e is a unit vector, e(z/,-) = 1. The BOOMERANG chan- 
nels consist of eight polarization-sensitive bolometers at 145 
GHz, and four spider-web bolometers at each of 245 GHz and 
345 GHz channels. We make use of data from all these detec- 
tors except two detectors that were known in pri or studies to 
be do minated by detector noise (245X and 345Z) dMasi et al. I 
120061) and two detectors with a significantly higher noise than 
the others (145Z2 and 345 Y), leaving us with seven detectors 
at 145, three at 245 and two at 345 GHz. 

We use the spectral response of each band as measured in 
lab with sub-percent accuracy. From these bands we derived 
the values of (0.49722, -0.21646, -1.01643) for s(v) at 145, 
245 and 345 GHz respectively. These bands provide ideal 
frequency coverage for an SZ study with channels in the SZ 
decrement, near the null, and the increment, respectively. The 
measured FWHM of the beams is 1 1 .5, 8.5, and 9. 1 arcmin for 
the 14 5, 245, and 345 GHz channel, respectively d Jones et al. I 
120061) . These values include a 2.4 arcmin pointing jitter. The 
beams window functions b'j a are in fact numerically derived 
from physical optics simulations, combined with a Gaussian 
pointing jitter. 

Similar to prior studies with BOOMERANG data, we pro- 
duced CMB temp erature anisotropy (T) maps using the Italian 
analysis pipeline dMasi et al. l 2006|) and the TT p ower spectra 
with the MASTER method dHivon et al. 112002b . To remove 
excess atmospheric noise, we filtered out time-ordered data 
in each of the detectors in frequency space below 200 mHz. 
This results in a damping of power at angular scales above 
1.2°. The effect of this filter, of the scanning strategy, of the 
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TABLE 1 



SZ Power Spectrum Estimates 





Bin 1 


Bin 2 


Bin 3 


^-range 


Ten /ICO 

zjU — 4jU 


4jU — /OU 


TOO 1 TOO 

/UU — IzUO 




Optimal wei 


ghts 




Vfl45GHz 


0.9323 


0.8514 


0.7289 


W245GHz 


0.4193 


0.3771 


0.3002 


H '345GHz 


-0.3515 


-0.2285 


-0.0292 


Raw SZ 


236 


164 


538 




Residuals 


(a) 




CMB 


53 


36 


70 


Instr. noise 


92 


12 


-95 


Galactic dust 


68 


82 


138 


FIRB 


44 


81 


195 


Radio sources 


3 


7 


58 


Total residual 


247 


202 


338 


SZ Band Power Uncertainties 


Instr. noise 


154 


116 


280 


Foregrounds 


37 


79 


145 


Beam 


3 


5 


44 


Calibration 


121 


77 


63 


Transfer func. 


2 


3 


11 


Cosmic & NG Var w 


7 


6 


4 


Final SZ Band Power 


-11 ± 199 


-38 ± 160 


200 ± 325 



NOTES. — The weights and Raw SZ designate the weight vectors for each 
multipole bin and the SZ power spectrum respectively with both as measured 
from data. Except in the case of weights w, the values are tabulated in units 
of fiK 2 for the SZ angular power spectrum / 2 C; /27r at the RJ end of the fre- 
quency spectrum. 

a: The residuals are the average spectra measured on our SZ-free simulations 
and represent our bias. The total residual is different from the sum of the par- 
tial residuals due to small (<10%) random correlation between components, 
b: The uncertainties are the dispersion measured with our simulations. The 
final SZ spectrum values are corrected for the noise and foreground bias with 
the dispersion error from simulations. 

c: Assuming the WMAP team's SZ power spectrum with erg = 0.95, the 
2(j upper limit we derived from all SZ data. The calculation for the Non- 
Gaussian (NG) covariance makes use of the same halo model as used for this 
power spectrum. 

sky coverage, and of the pixelization results in an effective 
window function. This window function is estimated by pro- 
jecting signal-only sky simulation into detector time-streams. 
To do that we used 10 sk y simulations of the SZ signal, ob- 
tained from I White! (120031) with different normalizations and 
initial conditions. These time-streams are then analyzed with 
our pipeline and the resulting angular power spectrum is com- 
pared to the power spectrum of the input sky signal to obtain 
the window function. The 10 different simulations result in 
slightly different window functions. The scatter in the win- 
dow function is included in the final error estimate of the SZ 
signal. 

3. SIMULATIONS 

Monte Carlo simulations were used in order to estimate the 
residual signal from correlated detector noise, primary CMB, 
and foregrounds that is detected as an SZ signal at the end. In 
fact, when we minimize the covariance in the data to extract 
the SZ, the optimal combination does not remove perfectly the 
other signals. This is particularly true with just three channels, 
and is expected to improve with more channels. We also use 
simulations to estimate final uncertainties in the SZ signal. 

We generate 200 time-stream simulations for each one of 
the 12 detectors with a combination of sky signal simulations 
that were projected into a time-line using the BOOMERANG 
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FIG. 2. — Angular power spectrum of SZ anisotropies at the RJ end of 
the frequency spectrum. We show the 68% confidence upper limits from the 
BOOMERANG data in hashed columns to the left and the reported upper 
limit at the 65% confidence level from the SZA experiment to the right. The 
data points from CBI, BIMA, and ACBAR experiments are also shown (see 
text for details). We correct the ACBAR point to RJ end of the frequencies. 
The lines show two theory predictions for SZ fluctuations: the curved region 
is a prediction based on results from numerical simulations with a^ z = 1 and 
the solid line is the same model as used by the WMAP team but scaled to 
o| z = 0.95, the 2-ct upper limit from all SZ data. The limit from BOOMER- 
ANG SZ data alone is o| z < 1.14 (95% c.l.). 

pointing, and a random noise realization with variance con- 
sistent with data. As the covariance matrix C,y of Eq. 1 is 
built using the cross-spectra between different detectors, the 
noise realizations include also correlated noise between dif- 
ferent detectors. For 14 5 GHz the noise co rrelation spectra 
are reported in Table 7 of Masi et al. I d2006l) and they are be- 
low 3% for the detectors used in the analysis. For the 245 and 
345 GHz we measured the noise correlation spectra removing 
the optimal signal map from the time-ordered data. The cor- 
relations are at most 6% for the detectors used in the analysis. 

In addition to 200 realizations of primary CMB, the sky sig- 
nal simulations are computed with Galactic dust, far-Infrared 
background (FIRB) sources, and radio point sources and 
combined with BOOMERANG bands passes. The Galac- 
tic d ust for the observed fie ld is described with model 8 
from lFinkbeiner et al. I (1 19991) . while FIRB and point sources 
were obtained from the Planck sky model (PSM: 1Eeach et al.~l 
120081: J. Delabrouille 2009, in preparation). 

The amplitude of each component in the data was estimated 
with a Monte Carlo Markov Chain for the three frequen- 
cies together. We compared the measured temperature angu- 
lar power spectrum with a linear combination of CMB, dust, 
FIRB, and point sources and fitted the amplitude coefficients 
p 1 of each component 7 such that Cj ot = ^ 7 /? 7 C ; 7 +Ni, where 
N/ is detector noise. The CMB amplitude of the predicted 
model is the same as we find in the data with /?cmb = 1.01 ± 
0.0 4. We find Pdust = 4.6 ± .8, larger than predicted by model 
8 of lFinkbeiner et al. I (1 19991) and p mliio = 0.76 ±0.07 for radio 
point sources, smaller than predicted bv lde Zotti et al. I (120051) 
counts used in PSM, but in agreement with iFriedman et al. I 
( 120091) . The amplitude of the FIRB remains unconstrained 
with < pfirb < 0.7 (lcr). 

To take the uncertainties in the amplitudes of foregrounds 
into account, we ran three sets of 200 simulations. In the first 
set we leave the level of foregrounds as the models predict 
(p a = 1 for each of CMB, dust, FIRB and radio sources). In 
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the second case we set pdust = 4.6, Phrb = 0, and p ra dio = 0.76. 
In the third case, we change /?firb = 0.7, while keeping the 
rest of the parameters as in the second case. We note that these 
simulations do not include the SZ signal as we are reconstruct- 
ing the SZ under the assumption of a zero signal. This does 
not bias the procedure since the optimal weights from Eq. 1 
are independent of the exact amplitude of the large-scale SZ 
effect. 

The simulated time-lines have been analyzed with the same 
pipeline as used for the data. This results in three sets of 
Monte-Carlo binned power spectra for the SZ Cf^ c , which 
should, in principle, be zero since the SZ signal is not in- 
cluded in the simulations. From the distribution of each set of 
Cf mc we derived ( 1 ) a bias in the SZ binned power spectrum 
from (Cf^o) and (2) the bin-to-bin covariance matrix due to 
statistical noise and sampling variance of the CMB, 



Cbb' - ((QmC (Qmc)) X ipb 



MC" 



' (^i' Z Mc))) ■ 



(3) 



The error bars on the angular power spectrum are given by 
ACf z = \JCbb- In Table [T] we quote the average value of the 
bias from the three sets in our residual amplitude and add the 
average of the dispersion of these residuals as an additional 
error (Table Q] foreground error). 

These foreground errors are combined with the beam er- 
rors, which are estimated again through Monte Carlo simu- 
lations. The calibration uncertainties of the time-lines were 
also included through Monte Carlo simulations. They are at 
most 2%, 8%, and 13% for the 145 GHz, 245 GHz, and 345 
GHz, respectively, leading to an error of 4%, 16% and 26% 
on the temperature angular power spectra at each of the three 
frequencies. All the uncertainties listed above are added in 
quadrature for the final SZ band-power uncertainty. 

4. SZ POWER SPECTRUM ESTIMATE 

In Table Q] we list the values we obtain for the three mul- 
tipole bins between I =250 and 1200. We also list the resid- 
ual level from each foreground component, including detector 
noise, and the error associated with various uncertainties as 
described above. As tabulated, the biggest contamination to 
SZ detection comes from instrumental noise at the largest an- 
gular scales, while the FIRB dominates the contamination at 
the smallest angular scales probed by the experiment. Ra- 
dio point sources generate negligible confusion, primarily 
because at these high frequencies radio sources produce a 
weaker background compared to the dusty galaxies making 
the FIRB. In accounting for foreground contamination in the 
SZ estimate, we have taken a conservative approach here al- 
lowing for all components. An aggressive approach with the 
assumption of no FIRB leads to a marginal detection of an SZ 
signal especially in the third bin. Though the amplitude of 
FIRB fluctuations at 350 GHz is uncertain and we have based 
our model on the PSM, we do not consider an SZ detection 
with a no FIRB assumption to be realistic. 

The binned SZ power spectrum limits at the RJ end of the 
frequency spectrum are shown in Figure [2] where we plot 
the 68% confidence level limit for three bins between multi- 
poles of 250 and 1200. In estimating the final SZ band power 
uncertainty, we also include the usual Gaussian cosmic vari- 
ance and the extra covariance from the non-Gaussian nature 



of the SZ power spectrum ( ICooravll2001h . This covariance is 
calculated assuming of 2 = 0.95 and making use of the same 
halo model as the one used by the W MAP team's SZ mode l 
and shown with a solid line in Fig. 2 dKomatsu et al. ||2~002|) . 
Within this model, we study the cosmological implications 
of our limit on the SZ fluctuations, using a MCMC package 
dLewis & Bridldl2002h to constrain the amplitude of fluctua- 
tions. For reference to numerical simulations, with a dashed 
line, we also show the average SZ sig nal and the sc atter from 
a set of ten simulations at a% = 1 from White ( 20Q3j). 

In addition to three bins shown in Figure [21 we also com- 
bine the estimation of BOOMERANG SZ power spectrum 
to a single broad bin of 250 < I < 1200. We find an up- 
per limit of 234 /xK 2 in 1(1+ 1)C//2tt at the 95% confidence 
level. Previous analytical calculations have shown that Cf z oc 
(af z ) 7 (Qbh) 2 dSeliaket al. 11200 lh . where we separate of z as- 
sociated with SZ from the primordial normalization a%. The 
amplitude constraint from BOOMERANG SZ data alone is 
of z < 1.14 at 95% confidence level. 

In Figure 12 we also compare our upper limits with re- 
sults on SZ fluctuat ions in th e literature, includin g CBI 
dSievers et al. I 120091). B IMA dDawson et al. I 120061) . and 
ACBAR dReichardt et al.1l2008l) . We scale the ACBAR value 
from 150 GHz to the RJ end of the spectrum for easy 
comparison with all other results. We also fit jointly the 
combined WMAP 5-y ear dKomatsu et al. 1 12008[). AC BAR 
dReichardt et al. 1 120081) . and CBI dSievers et al. I l2009h data 
together with SZ upper limits from BOOMERANG and SZA. 
We use the same analytical halo model with cr| z = 0.95 to 
include an extra uncertainty associated with non-Gaussian 
covariance in each of these measurements; these, however, 
make only a minor difference except in the case of BIMA 
where the smaller area surveyed increase the importance of 
non-Gaussianities. Marginalizing over all other cosmologi- 
cal parameters in the ACDM model, we find cr| z < 0.92 at 
95% confidence level (crf z < 0.71 at la). This SZ derived 
amplitude is fully consistent with WMA P 5-year result with 
erg = 0.81 ± 0.02 dKomatsu et al. 112008b . While it has been 
claimed in the past that the SZ derived erg is higher than the 
value derived from the CMB, we do not find this is the case 
with the BOOMERANG data. 

The next opportunities to perform a multi frequency anal- 
ysis similar to o urs will be with Planck and OLIMPO 
dMasi et aTl 120081) . Both these experiments include multi- 
ple bands at high frequencies where the SZ is positive. As 
we have found that only one channel above the SZ null fre- 
quency is not adequate to separate both CMB and FIRB from 
SZ fluctuations, with several high frequency channels, these 
upcoming CMB experiments should be able to measure and 
separate FIRB more accurately than we were able to with just 
one channel at 350 GHz. 
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